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Abstract — Q In this paper, collaborative use of relays to 
form a beamforming system and provide physical-layer security 
is investigated. In particular, amplify-and-forward (AF) relay 
beamforming designs under total and individual relay power 
constraints are studied with the goal of maximizing the secrecy 
rates when perfect channel state information (CSI) is available. 
In the AF scheme, not having analytical solutions for the optimal 
beamforming design under both total and individual power 
constraints, an iterative algorithm is proposed to numerically 
obtain the optimal beamforming structure and maximize the 
secrecy rates. Robust beamforming designs in the presence of 
imperfect CSI are investigated for decode-and-forward (DF) 
based relay beamforming, and optimization frameworks are 
provided. 

Index Terms: amplify-and-forward relaying, decode-and- 
forward relaying, physical-layer security, relay beamforming, 
robust beamforming. 

I. INTRODUCTION 

The broadcast nature of wireless transmissions allows for 
the signals to be received by all users within the commu- 
nication range, making wireless communications vulnerable 
to eavesdropping. The problem of secure transmission in 
the presence of an eavesdropper was first studied from an 
information-theoretic perspective in [ 1 1 where Wyner con- 
sidered a wiretap channel model. Wyner showed that secure 
communication is possible without sharing a secret key if 
the eavesdropper's channel is a degraded version of the 
main channel, and identified the rate-equivocation region and 
established the secrecy capacity of the degraded discrete 
memoryless wiretap channel. The secrecy capacity is defined 
as the maximum achievable rate from the transmitter to the 
legitimate receiver, which can be attained while keeping the 
eavesdropper completely ignorant of the transmitted messages. 
Later, Wyner's result was extended to the Gaussian channel in 
[|3l and recently to fading channels in [4] and [5|. In addition 
to the single antenna case, secrecy in multi-antenna models is 
addressed in (6l and Q. Regarding multiuser models, Liu et 
al. H| presented inner and outer bounds on secrecy capacity 
regions for broadcast and interference channels. The secrecy 
capacity of the multi-antenna broadcast channel is obtained 
in (9). Additionally, it is well known that even if they are 
equipped with single-antennas individually, users can cooper- 
ate to form a distributed multi-antenna system by performing 

1 This work was supported by the National Science Foundation under Grant 
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relaying. When channel side information (CSI) is exploited, 
relay nodes can collaboratively work similarly as in a MIMO 
system to build a virtual beam towards the receiver. Relay 
beamforming research has attracted much interest recently (see 
e.g., lfTOl -1 141 and references therein). Cooperative relaying 
under secrecy constraints was also recently studied in IfTSI - 
|18| . In [15 1, a decode-and-forward (DF) based cooperative 
protocol is considered, and a beamforming system is designed 
for secrecy capacity maximization or transmit power mini- 
mization. For amplify-and-forward (AF), suboptimal closed- 
form solutions that optimize bounds on secrecy capacity are 
proposed in |fl6l . However, in those studies, the analysis 
is conducted only under total relay power constraints and 
perfect CSI. In our recent work lfT8l . we studied the problem 
of DF beamforming under both total and individual power 
constraints. It is shown that the total power constraint leads 
to a closed-form solution. The design under individual relay 
power constraints is formulated as an optimization problem 
which is shown to be easily solved using two different 
approaches, namely semidefinite programming and second- 
order cone programming. In this paper, we study amplify- 
and-forward (AF) relay beamforming under both total and 
individual power constraints. We also extend our previous 
study by considering two robust beamforming design methods 
for DF relaying under imperfect CSI. 

II. Channel Model 

We consider a communication channel with a source S, a 
destination D, an eavesdropper E, and M relays {i? m }m=i 
as depicted in Figured] We assume that there is no direct 
link between S and D, and S and E. We also assume that 
relays work synchronously by multiplying the signals to be 
transmitted with complex weights {w rn } to produce a virtual 
beam point to destination and eavesdropper. We denote the 
channel fading coefficient between S and R m as g m G C, 
the fading coefficient between R m and D as h m € C, and 
the fading coefficient between R m and E as z m € C. In this 
model, the source S tries to transmit confidential messages to 
D with the help of the relays while keeping the eavesdropper 
E ignorant of the information. It's obvious that our channel is a 
two-hop relay network. In the first hop, the source S transmits 
x 8 to relays with power i?[|x s | 2 ] = P s . The received signal at 




Fig. 1. Channel Model 



R m is given by 



Vr,m fjrn^s T" ^rj 
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where r/ m is the background noise that has a Gaussian dis- 
tribution with zero mean and variance of N m . In first part of 
the paper, we focus on amplify-and-forward relaying. Robust 
design for DF will be discussed subsequently. 

In the AF scenario, the received signal at the m th relay 
R m is directly multiplied by l m w m without decoding, and 
forwarded to D. The relay output can be written as 



w, 



lrn{g II 



Vm)- 



The scaling factor, 
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is used to ensure £'[|ar T . jrrl | 2 ] = |w m | 2 . There are two kinds of 
power constraint for relays. First one is a total relay power 
constraint in the following form: ||w|| 2 = w^w < Pt where 
w = [wi, ...wm] T , Pt is the maximum total power for 
all relays. (-) T and (•)t denote the transpose and conjugate 
transpose, respectively, of a matrix or vector. In a multiuser 
network such as the relay system we study in this paper, 
it is practically more relevant to consider individual power 
constraints as wireless nodes generally operate under such 
limitations. Motivated by this, we can impose \w m \ 2 < •pra^rn 
or equivalently |w| 2 < p where | • | 2 denotes the element- wise 
norm-square operation and p is a column vector that contains 
the components {p m }- Vm is the maximum power for the m 
relay node. 

The received signals at the destination D and eavesdropper 
E are the superposition of the messages sent by the relays. 
These received signals are expressed, respectively, as 



M 



yd 



(g m Xs +Vm) +n , and 
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where no and n\ are the Gaussian background noise compo- 
nents at D and E, respectively, with zero mean and variance 



Nq. Now, it is easy to compute the received SNR at D and E 
as 



2~^m=l ^m^m^m^ml Ps 



E™=i \h m \ 2 l 2 m \w m \ 2 N n 
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and 



2~^m=l z mgm^m^m\ Ps 



Em=l \ z m\ 2 l 2 n \w m \ 2 Nr, 

The secrecy rate is now given by 

R s = I(x s ;y d ) - I(x s ;y e ) 

= log(l + T d ) - log(l + T e ) 
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I Em=l z mg m l m W m \ 2 P s + Em=l I 

\ 2 N m + N 
(10) 

where /(•; •) denotes the mutual information. In this paper, 
we address the joint optimization of {w m } with the aid 
perfect CSI and hence identify the optimum collaborative relay 
beamforming (CRB) direction that maximizes the secrecy rate 
in (QOI). 

III. Optimal Beamforming for AF case 
Let us define 

h s = [h* 1 g*l 1 ,...,h* M g* M l M } T , (ID 
h z = [ztglh, z* M g* M l M ] T , (12) 
D h = Diag(|/ii| 2 Z 2 Ai, \h M \ 2 l 2 M N M ), and (13) 
D z = ESag(|«i| 2 J?JVi, \z M \ 2 l 2 M N M ). (14) 

where superscript * denotes conjugate operation. Then, the 
received SNR at the destination and eavesdropper can be 
reformulated, respectively, as 



T d = 

r„ = 
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where tr(-) represent the trace of a matrix. It is obvious that we 
only have to maximize the term inside the logarithm function 
in ([Tol l. With these notations, we can write the objective 
function of the optimization problem as 



i + r d 
i + r e 

i + 
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wtD z w + A + P s wth z h z f w 

A + tr((D h + PshghgVwt) 

Ao + tr((D z +P s h z h z t )wwt) * A + tr(D h wwt) ' 

(19) 
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If we denote t\ — ^° +tr ( (D** +p° ) ww t j ' ^ 2 = l& obtained by solving following optimization problem: 



^°tw^! zWW -K , and define X = ww 1 " using the similar 
semidefinite programming method as described in |18|, we 
can express the optimization problem as 

max tii 2 

X,ti,t 2 

s.t tr(X(D z -t 2 D h )) >iV (t 2 -l) 

tr (X (p h + P.hghgt - h (D z + P s h z h z f ))) > JV (*i - 

and diag(X) < p, (and/or tr(X) < P T ) and X > 0. 

(20) 

where X y means that X is a symmetric positive semi- 
definite matrix. Since X by definition is a rank one matrix, 
finding the optimal weights is in general a nonconvex opti- 
mization problem. Thus, we above ignore the rank constraint, 
and hence employ semidefinite relaxation (SDR) lfl3l . If 
the matrix X opt obtained by solving the above optimization 
problem happens to be rank one, then its principal component 
will be the optimal solution to the original problem. 

Notice that this formulation is applied to both total relay 
power constraint and individual relay power constraint which 
are represented by tr(K) < Pt and diag(X) < p, respec- 
tively. When there is only total power constraint, we can easily 
compute the maximum values of t\ and f 2 separately since 
now we have Rayleigh quotient problems. 



max t\ 

s.t tr (X (D h + P s hgh g t - tx (D z + P s h z hJ 
and diag(X) < p, and X y 0, 

and 
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where A max (A, B) is the largest generalized eigenvalue of the 
matrix pair (A, B) 0. 

Similarly, maximum values of £ 2 under total power con- 
straint is 



t%u = max 



w f D z w + N 



tw<p T wtDhW + TVo 
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(24) 
(25) 
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When there are individual power constraints imposed on the 
relays, The maximum values t\^ u and ti,%,u El for t\ and 

2 For a Hermitian matrix A 6 C" x " and positive definite matrix B £ 
C" x ™, (A, 4>) is referred to as a generalized eigenvalue - eigenvector pair of 
(A, B) if (A, i>) satisfy Aip = XBtfi (T9). 

3 Subscripts i in t\ j u and t2.i,u are used to denote that these are the 
maximum values in the presence of individual power constraints. 



max t2 
x,t 2 



s.t tr (X (D z - i 2 D h )) > N (t 2 - 1) 
and diag(X) < p, and X h 0. 



) > N (h - 1) 
(27) 

(28) 



In fact, for any value of t%, the feasible set in (127t is convex. 
If, for any given ti, the convex feasibility problem 

find X 

s.t tr (X (D h + P.hghgt - ti (D z + P s h z h z ^ > N (ti - 1) 
and diag(X) < p, and X y 0, 

(29) 

is feasible, then we have ti t i :U > t\. Conversely, if the convex 
feasibility optimization problem ( |29l is not feasible, then we 
conclude ti : i M < ti. Therefore, we can check whether the 
optimal value ti t i tU of the quasiconvex optimization problem 
in d27| i is smaller than or greater than a given value ti by 
solving the convex feasibility problem d29t . If the convex 
feasibility problem ( |29l is feasible then we know > t. 

If the convex feasibility problem ( |29l is infeasible, then we 
know that ti t i lU < t\. Based on this observation, we can 
use a simple bisection algorithm to solve the quasiconvex 
optimization problem (l27j > by solving a convex feasibility 
problem ( 1291 at each step. We assume that the problem is 
feasible, and start with an interval [l,u] known to contain 
the optimal value ti t i )U . We then solve the convex feasibility 
problem at its midpoint ti = (I + u)/2 to determine whether 
the optimal value is larger or smaller than t. We update the 
interval accordingly to obtain a new interval. That is, if ti 
is feasible, then we set I — ti, otherwise, we choose u = t\ 
and solve the convex feasibility problem again. This procedure 
is repeated until the width of the interval is smaller than the 
given threshold. Then, we conclude that ti i u = I. 

Similarly, £2,1, u can be obtained with the same bisection 
algorithm by repeatedly solving the following feasibility prob- 
lems: 



find X 

s.t tr(X(D z -t 2 D h )) > JV (t 2 -l) 
and diag(X) < p, and X y 0. 



(30) 



To solve the convex feasibility problem, one can use the 
well-studied interior-point based methods as well. We use the 
well-developed interior point method based package SeDuMi 
[22 1, which produces a feasibility certificate if the problem is 
feasible, and its popular interface Yalmip |23|. 

Note that for both total and individual power constraints, 
the maximum values of ti and ti are obtained separately 
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above, and these values are in general attained by different 
X = ww^. Now, the following strategy can be used to obtain 
achievable secrecy rates. For those X values that correspond 
to t± t i tU and ti iU (i.e., the maximum t\ values under individual 
and total power constraints, respectively), we can compute 

the corresponding t 2 = ^ t ^Zt) and denote them 
as t2,i,i and t 2 j for individual and total power constraints, 
respectively. Then, log(tx t i tU t2,i,i) and log(ti tU t 2t i) will serve 
as our amplify-and-forward achievable rates for individual and 
total power constraints, respectively. With the achievable rates, 
we propose the following algorithm to iteratively search over 
t\ and t 2 to get the optimal t\ i0 and t 2 ^ that maximize the 
product t\t 2 by checking the following feasibility problem. 

find X y 

s.t tr (X (D, - t a D h )) > N Q {t 2 - 1) 

tr (x (D h + P.hghgt - h (D, + P s h z h z W > N (h - 1) 

and tr(X) < Pt if there is total power constraint, 
or diag(X) < p if there is individual power constraint. 



IV. Robust Beamforming Design for DF case 

In the second hop, we can also employ decode and forward 
transmission scheme. In this scheme, each relay R m first 
decodes the message x s and normalizes it as x' s = x s /^/P^. 
Subsequently, the normalized message is multiplied by the 
weight factor w m to generate the transmitted signal x r = 
w m x' s . The output power of each relay R m is given by 

|2l _ 77N„,, J |2l _ |„„ |2 



A. Proposed Algorithm 
Define the resolution At 



N 



or At 



N 



(31) 



for some 



large N for total and individual power constraints, respectively. 

1) Initialize ii = ti iU , t 2 . Q — t 2 j when total power 
constraint is imposed, and ii = t± t i tU , t 2t0 — t 2 ,i.i 
when individual power constraint is imposed. Initialize 
the iteration index i = N. 

2) Set t\ = iAt. If tit 2tU < t\^ t 2y0 (total power constraint) 
or tit 2t i, M < ti, t 2 , (individual power constraint), then 
go to Step (3). Otherwise, 



a) Let t 2 



ti,o*a„ 



Check the feasibility problem 



OTb . If it is infeasible, go to step (3). If it is 
feasible, use the bisection algorithm in OTT ) with 
t% to get the maximum possible values of t 2 and 
denote this maximum as t 2 . m . The initial interval 
in the above bisection algorithm can be chosen as 



rt ll0 t 2 „ 



[t 1:0 t 2 , 



, t 2 ,i, u ] depending on the 



i = i — 1. Go back 



power constraints, 
b) Update £i j0 = t\, t 2:0 = t 2 , 
to step (2). 

3) Solve the following convex problem to get the optimal 



mm 
x 



tr(X) 

s.t tr (X (D z - i 2 , D h )) > N {t 2 . o - 1) 

tr (X (D h + Pshghgt - t lt0 (D z + P s h z h,t ^ 

> iVo(ti, - 1) 
X h and 

tr(X) < Pt if there is total power constraint, 

diag(X) < p if there is individual power constraint. 

(32) 



E[\x r \ 2 ]=E[\w m x'X] = \WmY 



(33) 



The received signals at the destination D and eavesdropper 
E are the superpositions of the signals transmitted from the 
relays. These signals can be expressed, respectively, as 

M 

y d = ^2 h m w m x' s + n = h f wa;^ + n , and (34) 

m— 1 
M 



I), = ^ Z rnW m x' s +n 1 = z'w^ + U x . 



(35) 



Additionally, we have defined h = [/if, ....h* M ] T , z = 
[z* , ....z^y, The metrics of interest are the received SNR 
levels at D and E, which are given by 



r l Em=l h m w m \ 2 I Eggl z mW m \ 2 

Td = N andre = W a ' (36) 

It has been proved that the secrecy rate R s over the channel 
between the relays and destination is 



R s = I{x s ;y d ) - I(x s ;y e ) 

= bg(i + r d ) - io g (i + r e ) 



log 



T M h w |2 " 

/ ; r — i * brawn 



.^0+ \ Y.m = l Z rnW. 



(37) 
(38) 

(39) 



By using the semidefinite relaxation (SDR) approach to ap- 
proximate the problem as a convex semidefinite programming 
(SDP) problem, the beamforming design for decode-and- 
forward under perfect CSI with individual power constraint 
can be formed as the following optimization problem [18]: 



max t 
x,t 



s.t tr(X(hh t - izz+)) > N (t - 1), 
and diag(X) < p, and X ^ 0. 



(40) 



which can be solved efficiently by interior point methods with 
a bisection algorithm fl8l . 

Systems robust against channel mismatches can be ob- 
tained by two approaches. In most of the robust beamforming 
methods, the perturbation is modeled as a deterministic one 
with bounded norm which leads to a worst case optimization. 
The other approach applied to the case in which the CSI 
error is unbounded is the statistical approach which provides 
the robustness in the form of confidence level measured by 
probability. 

We define H = hlr and Z = zz^ as the channel estimators, 
and H = H — H and Z = Z — Z as the estimation errors. 
First, consider the worst case optimization. In the worst case 
assumption, H and Z are bounded in their Frobenius norm as 
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||H|| < eh, ||Z|| < ez, where eh, £z are assumed to be upper 
bounds on the channel uncertainty. Based on the result of l20ll . 
the robust counterpart of the previously discussed optimization 
problem can be written as 

max t 
x,t 

s.t tr(X((H-e ff I)-i(Z + e z I)) >N (t-l), < 41 > 
and diag(X) < p, and X y 0. 

Note that the total power constraint tr(X) < Pt can be added 
into the formulation or substituted for the individual power 
constraint in (|4TT >. This problem can be solved the same way 
as discussed in fl8l . 

However, the worst-case approach requires the norms to be 
bounded, which is usually not satisfied in practice. Also, this 
approach is too pessimistic since the probability of the worst- 
case may be extremely low. Hence, statistical approach is a 
good alternative in certain scenarios. In our case, we require 
the probability of the non-outage for secrecy transmission is 
greater than the predefined threshold e by imposing 



Pr [ ^o+tr((H + H)X) > ^ 



N + tr((Z + Z)X) 



= Pr (tr (x(H + H - t(Z + Z)) >N (t- 1)) ) > e. 



(42) 



Now, the optimization problem under imperfect CSI can be 
expressed as 

max t 
x,t 

s.t Pr (tr (x(H +H - t(Z + Z)) > N (t -1))) > e, 
and diag(X) < p (or tr(X) < P T ), and X > 0. 

(43) 

If relays are under individual power constraints, we use 
diag(K) < p. Otherwise, for the case of total power con- 
straint, we use ir(X) < Pr- We can also impose both 
constraints in the optimization. For simplicity of the analysis 
we assume that the components of the Hermitian channel es- 
timation error matrices H and Z are independent, zero-mean, 
circularly symmetric, complex Gaussian random variables with 
variances <r|j and a 2 ^. Now, we can rearrange the probability 
in the constraint as 



Pr (tr ((H - tZ + H - tZ)X\ >(t- 1)N ] 



(44) 



Let us define y = tr ^(H — tZ + H— tZ)Xj . For given 

X, H, and Z, we know from the results of fl2~Tl that y 
is a Gaussian distributed random variable with mean \i = 



Then, the non-outage probability can be written as 

i ( (y-^) 2 



Pr(y >(t- 1)N ) 



(t-l)JVo 



27T(T, 



■ exp 



2al 



(45) 



1-lerf Wo-* ) >£ 
2 2 I y/2a v J ~ ' 



or equivalently as, 

(t - 1)N - ji 
\/2a y 



< err 1 (-2e + 1) 



(46) 



(47) 



Note that e should be close to one for good performance. 
Thus, both — 2e + 1 and ( t ~ 1 ) jv Q~^ gho^d ^ e negative valued. 

V2(7y ° 

Note further that we have tr(XX^) = ||X|| 2 , and hence cr y = 

an be written a; 

(t - 1)N - n 



a 2 - + t 2 cr|||X||. Then, this constraint can be written as 



X < 



2(4+tV|)erf- 1 (-2 £ + l) 



(48) 



As a result, the optimization problem becomes 

(t - 1)N - M 



max t 



s.t 11X11 < 



/2(4 + i 2 ( r|)err 1 (-2£ + l)' 
and diag(X) < p(or tr(X) < P T ), and X h 0. 



(49) 



tr ( (H — £Z)X) and variance a 2 



(4+tV|)tr(XXt) 



Using the same bisection search, we can solve this optimiza- 
tion numerically. 

V. Numerical Results 

We assume that {g m }, {h m }, {z m } are complex, circularly 
symmetric Gaussian random variables with zero mean and 
variances cr 2 , cr 2 , and tr 2 respectively. Moreover, each figure 
is plotted for fixed realizations of the Gaussian channel coeffi- 
cients. Hence, the secrecy rates in the plots are instantaneous 
secrecy rates 

In Fig. |2] we plot the secrecy rate for amplify-and-forward 
collaborative relay beamforming system for both individual 
and total power constraints. We also provide the result of 
suboptimal achievable secrecy rate for comparison. The fixed 
parameters are a g = 10, ah = 2, <r z = 2, jV m = l,iVo = 1 an d 
M = 10. Since the AF secrecy rates depend on both the source 
and relay powers, the rate curves are plotted as a function of 
Pt/P s - We assume that the relays have equal powers in the 
case in which individual power constraints are imposed, i.e., 
Pi = Pt/M. It is immediately seen from the figure that the 
achievable rates for both total and individual power constraints 
are very close to the corresponding optimal ones. Thus, the 
achievable beamforming scheme is a good alternative in the 
amplify-and-forward relaying case due to the fact that it has 
much less computational burden. Moreover, we interestingly 
observe that imposing individual relay power constraints leads 
to only small losses in the secrecy rates with respect to the 
case in which we have total relay power constraints. 
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Fig. 3. DF second secrecy rate vs. Pt under different e. 



In Fig. [3] we plot the maximum second hop secrecy rate of 
decode-and-forward that we can achieve for different power 
Pt and non-outage probability e values. In this simulation, 
we fix M = 5. h and z are randomly picked from Rayleigh 
fading with <j- h = 1 and <Jz = 2, and we assume that estimation 
errors are inversely proportional to Pt- More specifically, in 
our simulation, we have <j 2 ^ = 0.1/ Pt and cr| = 0.2/ Pt- We 
also assume the relays are operating under equal individual 
power constraints, i.e., pi — It is immediately observed in 
Fig. [3] that smaller rates are supported under higher non-outage 
probability requirements. In particular, this figure illustrates 
that our formulation and the proposed optimization framework 
can be used to determine how much secrecy rate can be 
supported at what percentage of the time. For instance, at 
Pt = 100, we see that approximately 7 bits/symbol secrecy 
rate can be attained 70 percent of the time (i.e., e = 0.7) 
while supported secrecy rate drops to about 5.8 bits/symbol 
when e = 0.95. 
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